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CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] This application claims the benefit under 35 U.S.C. Section 1 19(e) of the 
following co-pending and commonly- assigned U.S. provisional patent application(s), 
which is/are incorporated by reference herein: 
5 [0002] Provisional Application Serial No. 60/47 1 ,803, filed on May 20, 2003, by 
Mladen Barbie, entitled "Method for High Resolution Two-Dimensional 
Tomographic Magnetic Resonance Imaging Using Magnetic Probes," attorneys' 
docket number 176.19-US-P1 (CIT-3891-P). 

10 STATEMENT REGARDING FEDERALLY SPONSORED 

RESEARCH AND DEVELOPMENT 
[0003] The invention was made with Government support under Grant No. DAAD 
19-00-1-0392 awarded by the Army Research Office and Grant No. F49620-01-1- 
0497 awarded by the Air Force Office for Scientific Research. The Government has 

1 5 certain rights in this invention. 

BACKGROUND OF THE INVENTION 
1 . Field of the Invention. 

[0004] The present invention relates generally to magnetic resonance imaging, and 
20 in particular, to a method, apparatus, and article of manufacture for high resolution 
two-dimensional tomographic magnetic resonance imaging using magnetic probes. 
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2. Description of the Related Art. 

[0005] (Note: This application references a number of different publications as 
indicated throughout the specification by one or more reference numbers within 
brackets, e.g., [x]. A list of these different publications ordered according to these 
5 reference numbers can be found below in the section entitled "References." Each of 
these publications is incorporated by reference herein.) 

[0006] Magnetic resonance imaging [1][2] (MRI) has advanced at a rapid pace since 
initial proposals and demonstrations in 1973 [3] [4], with applications in medical 
imaging attracting the most attention [5] [6]. Advances in magnetic resonance 

10 microscopy? have also been significant [8], recently reaching two-dimensional (2D) 
imaging resolution of 1 ^m [9]. Improvements in conventional inductive methods of 
magnetic resonance detection [10][1 1] and application of imaging gradients have 
generally been used to achieve such advances. Although the possibilities of 
angstrom-scale resolution were pondered in the early days of MRI [4] [12], the 

15 ultimate goal of achieving atomic resolution has remained elusive. 

[0007] In 1991, an alternative method of applying the imaging gradients and 
detecting magnetic resonance, magnetic resonance force microscopy (MRFM)[13] 
was proposed, with the ultimate goal of single spin sensitivity and three-dimensional 
(3D) imaging capability. The technique relies on the atomic scale imaging gradients 

20 from the microscopic magnetic particle mounted on a micromachined mechanical 

cantilever for the appropriate detection sensitivity required for 3D single spin imaging 
[14]. Successful MRFM demonstrations were reported for the cases of electron spin 
[15] nuclear spin [16], and ferromagnetic [17] resonance systems. MRFM research 
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has benefited from the low temperature implementations of the instrument [18], and 
rapid advances in the fabrication techniques for incorporating smaller magnetic 
particles [19] [20], and more sensitive mechanical resonators [21]. However, reported 
MRFM imaging resolution of -1 jam [22] [23] remains at the level of inductive 
5 detection in conventional MRL 

[0008] A complementary atomic resolution magnetic resonance imaging method 
[24] may be used to significantly relax the challenging technical requirements of 
single spin detection by imaging discrete ordered crystal lattice planes where many 
spins coherently contribute to the magnetic resonance signal. Such an approach 

10 closely resembles the initial magnetic resonance imaging proposal [4] [12] in which 
linear magnetic field gradients are used to selectively excite magnetic resonance in 
different atomic lattice planes and produce "diffraction"-like effects. However, the 
approach differs from the initial magnetic resonance imaging proposal by introducing 
nonlinear magnetic fields and field gradients from a ferromagnetic sphere to achieve 

15 atomic resolution magnetic resonance diffraction. By using various sample-detector 
coupling mechanisms [25], it may be shown that the realization of the long desired 
atomic resolution magnetic resonance diffraction of crystals may be within reach of 
available experimental techniques. 

[0009] However, the prior art lacks the ability to conduct two-dimensional magnetic 
20 resonance microscopy of noncrystalline structures at an atomic resolution. 

SUMMARY OF THE INVENTION 
[0010] The invention introduces the concept of computerized tomographic 
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microscopy in magnetic resonance imaging using the magnetic fields and field 
gradients from a ferromagnetic probe. A configuration may be used wherein a sample 
(two-dimensional or three-dimensional) is under the influence of a large static 
polarizing field, a small perpendicular radio-frequency field, and a magnetic field 
5 from a ferromagnetic sphere. 

[0011] Despite the nonuniform and nonlinear nature of the fields from a 
microscopic magnetic sphere, the concepts of computerized tomography can be 
applied to obtain proper image reconstruction from the original spectral data by 
sequentially varying the relative sample-sphere angular orientation. Accordingly, 
10 atomic resolution magnetic resonance imaging of discrete periodic crystal lattice 
planes using ferromagnetic probes may be extended to two-dimensional imaging of 
noncrystalline samples with resolution ranging from micrometer to angstrom scales. 

BRIEF DESCRIPTION OF THE DRAWINGS 
1 5 [0012] Referring now to the drawings in which like reference numbers represent 
corresponding parts throughout: 

[0013] FIG. 1 shows the contours of constant values of the z component of a 
magnetic field from a ferromagnetic sphere in accordance with one or more 
embodiments of the invention; 
20 [0014] FIGS. 2(a), 2(b), and 2(c) reproduce the magnetic resonance spectra for 
crystalline structures next to a 100 nm diameter cobalt ferromagnetic sphere for the 
three cases: (a) semiinfinite crystal, (b) thin film 100-unit-cells in thickness, and (c) 
100-by-lOO-by-lOO atoms crystallite in accordance with one or more embodiments of 
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the invention; 

[0015] FIGS. 3(a), 3(b), and 3(c) show three representation plots for crystal lattice 

spin sites in accordance with one or more embodiments of the invention; 

[0016] FIGS. 4(a), 4(b), and 4(c) illustrate the positioning of a two-dimensional 
5 sample with respect to the ferromagnetic sphere and the sphere's magnetic field 

contours in accordance with one or more embodiments of the invention; 

[0017] FIG. 5 shows the schematic representation of the configuration and 

accompanying parameters that may be used in an image reconstruction process in 

accordance with one or more embodiments of the invention; 
10 [0018] FIGS. 6(a)-6(d) illustrate a two-dimensional phantom simulated at two 

different angles (FIGS. 6(a) and 6(c)) and accompanying magnetic resonance spectra 

(FIGS. 6(b) and 6(d)) in accordance with one or more embodiments of the invention; 

[0019] FIG. 7(a) illustrates a projection at 45° superimposed by an ideal projection 

in accordance with one or more embodiments of the invention; 
15 [0020] FIG. 7(b) illustrates a reasonable reconstruction of a sample spin density 

without linearization and de-convolution, in accordance with one or more 

embodiments of the invention; 

[0021] FIG. 8(a) illustrates a linearized projection at 45° superimposed by an ideal 
linearized and de-convolved projection in accordance with one or more embodiments 
20 of the invention; 

[0022] FIG. 8(b) illustrates a reasonable reconstruction with non-linear distortions 

eliminated in accordance with one or more embodiments of the invention; 

[0023] FIG. 9(a) illustrates a linearized and de-convolved projection in accordance 




with one or more embodiments of the invention; 

[0024] FIG. 9(b) illustrates an image reconstruction in accordance with one or more 
embodiments of the invention; 

[0025] FIGs 10(a)-10(c) show the computed images of three successively improved 
5 reconstructions, respectively, where 90 simulated projections at 2° apart are used in 
accordance with one or more embodiments of the invention; 
[0026] FIG. 1 1 illustrates several possible sample-sphere orientations and the 
respective sample rotation axis required for computerized tomographic image 
reconstruction in accordance with one or more embodiments of the invention; 
10 [0027] FIG. 12 illustrates an example device wherein a coil is positioned near a 

ferromagnetic sphere on the cleaved end of the optical fiber in accordance with one or 
more embodiments of the invention; 

[0028] FIG. 13 illustrates a sample device wherein micro-coil structures are wound 
on capillary tubes filled with ferromagnetic material in accordance with one or more 

1 5 embodiments of the invention; 

[0029] FIGS. 14(a) and 14(c) illustrate electron microscope images of a multi- 
functional nanostructure/nanowire resonator that may be used for high sensitivity 
detection in accordance with one or more embodiments of the invention; 
[0030] FIG. 14(b) illustrates an optical microscopic image of the nanowire resonator 

20 of FIG. 14(a) in accordance with one or more embodiments of the invention; 

[0031] FIG. 14(d) illustrates the method of suspending the nanowire resonator of 
FIGS. 14(a), (b), and (c) in accordance with one or more embodiments of the 
invention; 
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[0032] FIG. 15 illustrates an interferometric optical detection set-up that could be 
used for the measurement of a nano-mechanical system in accordance with one or 
more embodiments of the invention; and 

[0033] FIG. 16 is a flow chart illustrating the logical flow for conducting magnetic 
resonance tomographic microscopy in accordance with one or more embodiments of 
the invention 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0034] In the following description, reference is made to the accompanying 
drawings which form a part hereof, and which is shown, by way of illustration, 
several embodiments of the present invention. It is understood that other 
5 embodiments may be utilized and structural changes may be made without departing 
from the scope of the present invention. 

Overview 

[0035] In the present invention, a method extends the use of imaging gradient 
10 magnetic fields from ferromagnetic probes to magnetic resonance microscopy of non- 
crystalline disordered samples. Despite the non-uniform and nonlinear nature of the 
fields from a microscopic magnetic sphere, the concepts of computerized tomography 
can be applied to obtain proper image reconstruction from the original spectral data 
by sequentially varying the relative sample-sphere angular orientation. 
1 5 [0036] The description below first reviews the concept of atomic resolution 
magnetic resonance diffraction of ordered crystals using the fields from a 
ferromagnetic sphere in order to set the stage for extending the technique to 2D 
magnetic resonance tomographic microscopy of non-crystalline samples. 

20 Magnetic Resonance Imaging of Discrete Crystal Lattice Planes 

[0037] The previously described model of atomic resolution magnetic resonance 
imaging of crystal lattice planes [24] [25] relies on placing a ferromagnetic sphere in 
proximity of the surface of a crystal. A cobalt sphere, ro=50 nm in radius, has a 
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magnetization per unit volume of 1500 emu/cm 3 . The simple cubic lattice crystal is 
assumed to have a unit cell size of ao=3 A. A large dc magnetic field Bo is applied 
parallel to the sample surface in the z direction, polarizing the spins of the atomic 
lattice as well as saturating the magnetization of the ferromagnetic sphere. A small 
5 radio frequency field B\ is applied perpendicular to the large polarizing dc magnetic 
field Bo. The magnetic field from the ferromagnetic sphere at point r in the sample 
has the following azimuthally symmetric dipolar form: 

B(r) = ' , [Equation ( 1 )] 

w 

where n is the unit vector that points from the center of the ferromagnetic sphere to 
10 the crystal site location and m is the magnetic moment vector of the sphere. 

[0038] Since the external dc polarizing magnetic field Bo is considered to be much 
larger than the field from the ferromagnetic sphere, only the z component of the 
magnetic field from the ferromagnetic sphere is included when considering the 
resonant spins of the atomic lattice [7][24]-[26]: 

M 

15 = — f(3cos 2 0-l), [Equation (2)] 

M 

where #is the angle between the z axis and the distance vector r, and Mo is the 

magnitude of the saturation magnetic moment of the ferromagnetic sphere. 

[0039] FIG. 1 shows the contours of constant values of the z component of the 

magnetic field from the sphere, Bz , that have the azimuthally symmetric form around 

20 the z axis. As illustrated, the spins of the sample laying on the same contour have the 

same magnetic resonance frequency. For a lOOnm diameter cobalt sphere, the 

magnetic field gradients are sufficiently large that the discrete nature of the spins of 
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the sample may need to be taken into account. 

[0040] In the absence of the ferromagnetic sphere, the discrete nuclear spin sites in 
the crystal may experience the same externally applied field Bo and therefore meet the 
magnetic resonance condition at the same magnetic resonance frequency co/?. 
5 However, close to the ferromagnetic sphere, a large magnetic field gradient penetrates 
into the crystal, and only certain spin sites of the lattice satisfy the correct magnetic 
resonance conditions at any given magnetic field and frequency values. 
[0041] In the initial model [24][25], a numerical summation was computed to 
construct a histogram of the number of resonant spin sites in the sample within a 1 G 

10 wide shell of constant Z?z. This value of the bin width was selected since the linewidth 
broadening in solids is on the order of 1 G [26]. Distinct spectral peaks were 
discovered in the number of resonant spin sites with respect to the applied magnetic 
field in the negative value range. FIGS. 2(a), 2(b), and 2(c) reproduce the magnetic 
resonance spectra for crystalline structures next to a 100 nm diameter cobalt 

15 ferromagnetic sphere, between the field range of ifo—lOOOG and 5o=-500G for the 
three cases: (a) semiinfmite crystal, (b) thin film 100-unit-cells in thickness, and (c) 
100-by-lOO-by-lOO atoms crystallite, with the insets of the figures indicating the 
sample-sphere relative positions. A reduction of the sample size results in increased 
spectral peak contrast. 

20 [0042] The appearance of spectral peaks was explained using the 3D plots of the 
resonant spins under the influence of the polarizing magnetic field Bo and the 
magnetic field from a ferromagnetic sphere [24][25]. FIGS. 3(a), 3(b), and 3(c) show 
three such representation plots for crystal lattice spin sites that are in resonance at the 
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magnetic field value of i?o-625 G, the location of one of the sharp resonant peaks in 
the spectra of FIGS. 2(a), 2(b), and 2(c). A black dot is placed for every lattice site 
for which the field at that site is within 1 G of the resonant magnetic field. Only the 
positive values for the y axis are plotted for clarity. Dark bands of spins from the 
5 crystal lattice planes perpendicular to the magnetic field direction are responsible for 
the sharp spectral peaks. For semi-infinite crystal in FIG. 2(a), there is a large 
background signal that is diminished for the thin crystalline film in FIG. 2(b) since 
there are no resonant spins beyond the 100 th cell in the x-direction. The background 
signal is further minimized for the case of a small crystallite in FIG. 2(c). The 
10 discrete nature of the spins in the crystal lattice is required for the observation of 
magnetic resonance spectral peaks. 

[0043] At this magnetic field value (e.g., at -625 G), sections of the shell of 
constant Bz perpendicular to the crystal surface intersect the crystal lattice so that a 
large number of spin sites from the lattice plane at the top and bottom sections of the 
15 resonant shell satisfy the resonance condition. The bands of the resonant atoms from 
the lattice planes are clearly visible, and the resonant ring bands like these at the 
distinct magnetic field values are responsible for the sharp peaks in the magnetic 
resonance spectra of FIGS. 2(a), (b), and (c). 

[0044] It is emphasized that the appearance of the sharp magnetic resonance 
20 spectral peaks may only be possible by incorporating into the model the discrete 
nature of the atomic lattice sites, since magnetic resonance of a continuous medium 
would result in the monotonic spectrum. It is also noted that this magnetic resonance 
imaging technique specifically requires Bo to be significantly larger than the field 
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from the ferromagnetic sphere, say -10 (T) external magnetic field versus 0.1-1 (T) 
field from the ferromagnetic sphere. If the external magnetic field is comparable in 
magnitude to the magnetic field from the ferromagnetic sphere, contours of constant 
magnetic field magnitude may need to be considered. A review of magnetic fields in 
5 this regime [27] reveals that the contours of constant magnitude of magnetic field may 
not have the appropriate azimuthally symmetric form for the technique to be viable. 

Tomographic Magnetic Resonance Microscopy 

[0045] Following the description of imaging of crystal lattice planes using 
10 ferromagnetic probes, the main question remains as to whether the technique can be 
extended to the imaging of noncrystalline samples. Inspection of the 3D images of 
FIGS. 3(a), 3(b), and 3(c), showing the resonant spin sites at a particular value of the 
magnetic field for crystalline samples with different dimensions, reveals particular 
features that open the possibility of extending the technique to tomographic imaging 
15 of noncrystalline samples. 

[0046] It is apparent from the spectrum of FIG. 2(a) and the 3D image of FIG. 3(a) 
that the sharp spectral peaks for the semi-infinite crystal come from the very narrow 
regions of the sample, while there is a large background signal from other resonant 
spin sites that are intersected by the 1-G-thick shell of constant Bz . For the thin film 
20 with 100-unit-cell thickness of FIG. 2(b), the spectral contrast is significantly 
increased. The resonant peaks occur at the same location as for the semi-infinite 
crystal of FIG. 2(a), but they lack the large background signal because there are no 
atoms intersected by the resonant shell beyond the 100th unit cell along the x 




direction, as shown in the FIG. 3(b). 

[0047] There is still a background signal in FIG. 2(b) from the spins that are not part 
of the crystal lattice planes of interest, but their effect on the spectrum contrast is 
significantly reduced. The background signal from the spins away from the crystal 
5 lattice planes of interest are further minimized for the case of the small crystallite 
shown in FIGs. 2(c) and 3(c). 

[0048] Since the spectral contrast may be increased by the reduction of the sample 
size, the idea of tomographic magnetic resonance imaging of two-dimensional 
samples may be introduced. The two dimensional samples have sizes of ~ 1/10 the 
10 size of the ferromagnetic sphere and may be positioned as shown in FIGS. 4(a), 4(b), 
and 4(c). 

[0049] The sample can represent either a small molecule or protein -lO nm in size 
underneath a 100 nm diameter ferromagnetic sphere, or a biological cell --10 jam in 
size under a sphere 100 j^m in diameter. The use of the sphere model and these 

15 sphere dimensions are reasonable, as ferromagnetic spheres in this size range have 
been successfully fabricated [28]— [31] and have already been integrated as probes on 
cantilever structures for unrelated applications [32]-[34]. Modeling of the technique 
that includes sphere imperfections (that would undoubtedly introduce imaging 
aberrations) may also be conducted. 

20 [0050] FIG. 4(a) illustrates the positioning of the ferromagnetic sphere in relation to 
the sample. By magnifying the side view of the model arrangement, as shown in FIG. 
4(b), it is observed that the sample is intersected by approximately perpendicular lines 
of constant Bz along the x direction. A second view is shown in FIG. 4(c) where the 
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front view of the configuration of FIG. 4(a) is shown. By taking into consideration 
the large radius of curvature of the azimuthally symmetric intersecting contours of 
constant Bz (shown in FIG. 4(a)), it is apparent that the sample is intersected by 
approximately parallel lines of constant Bz along the y direction. Therefore, for an 
5 approximately flat two-dimensional sample with lateral dimensions smaller than the 
ferromagnetic sphere size, the magnetic resonance spectra will be the one- 
dimensional projections of the spin density along the z axis. 

[0051] However, as FIG. 4(c) shows, the lines of constant Bz are not equally spaced 
along the z axis, since the magnetic field is not linear along the z axis. Nevertheless, 
10 since the magnetic field is well defined in the sample region (see Equation (2) above), 
this nonlinearity can be compensated for in the computerized tomographic image 
reconstruction. 

[0052] The realization that the magnetic resonance spectrum of the sample in the 
configuration shown in FIG. 4 is a one-dimensional projection of the spin density 
15 along the z axis leads to the possibility of performing a high-resolution computerized 
tomographic reconstruction imaging. 

[0053] Imaging of samples through projections [35] has been an important concept 
ever since the discovery of x-rays [36] and has been used in gravitational theory [37] 
and radio astronomy [38] before becoming widespread through computerized 
20 tomography [39] in x-ray [40], electron [41], and optical [42] imaging, among others. 
Computerized tomographic image reconstruction algorithms are well known [43]- 
[45]. One or more embodiments of the present invention may apply similar principles 
to the case of magnetic resonance tomographic microscopy using ferromagnetic 
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spheres. However, it should also be noted that the image reconstruction from 
projections in magnetic resonance dates back to the first report of spin distribution 
imaging [3], and is still performed in the technique of stray field magnetic resonance 
imaging [46] where constant magnetic field gradients, on the order of 60 T/m, from 
5 superconducting magnets are used. 

[0054] One or more embodiments of the invention provide for the use of ultrahigh 
gradient fields (from microscopic ferromagnetic probes, as used in magnetic 
resonance force microscopy [14] [-5 X 10 6 T/m for a 100 nm diameter cobalt 
sphere]), in tomographic imaging of non-crystalline samples with resolution reaching 

10 angstrom levels. 

[0055] FIG. 5 shows the schematic representation of the configuration with the 
parameters used in the image reconstruction process in accordance with one or more 
embodiments of the invention. In the conventional computerized tomography [41], 
uniformly separated parallel rays are used to obtain an image projection along an axis, 

15 and the one-dimensional Radon transform of the sample density function p(y^) is 
formed 

(t,Q)imeP(y> z )d s • [Equation (3)] 

[0056] By obtaining a multiple of one-dimensional Radon transforms (Equation (3)) 
at different angles fa image reconstruction is performed by the Fourier transform 
20 filtered back projection algorithm for parallel projections [7]: 

= r{C[C P ^^^]l*^ 2 *^}^ [Eq^on (4)] 
[0057] This reconstruction process involves calculation of the Fourier transform of 
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the radon transform (innermost bracketed term), multiplication by a ramp function \k\ 
in conjugate space followed by an inverse transformation (outer bracketed term), and 
finally integration over all angles for the completion of the image reconstruction 
(outermost integration term) [7]. 
5 [0058] In one or more embodiments of the invention, a two-dimensional phantom is 
simulated with equally separated 9-by-9 array of spins in the y - z plane as illustrated 
in FIGs. 6(a)-6(d). It may be initially assumed that the points are individual proton 
nuclear spins separated by 5 A, each having a Lorentzian resonant response with a 10 
G linewidth. In the Cartesian coordinate system, the z component of the magnetic 
10 field from a ferromagnetic sphere [Equation (2)] has the form 

B^y,z)^ M f Z \ x2 ;^ ] [Equation^)] 
(x + y +z ) 

[0059] The center of the phantom array is placed next to the sphere as shown in 
FIG. 4, so that the contours of constant Bz of Equation (5) are approximately 
perpendicular to the phantom along the x axis, and parallel to the phantom along they 
15 axis. This location is determined by the condition that 

^=A = 0 [Equation (6)] 

ox oy 

and it was previously shown that at this location [25] 

x 2 + y 2 =4z 2 [Equation (7)] 
[0060] Therefore, when the sample is centered on they axis (yo=0), as shown in 
20 FIG. 4(c), the phantom is positioned along the z axis as shown in FIGs. 4(a)-4(b) at 
the 20 value of 

r 0 2 = 4zl . [Equation (8)] 
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[0061] FIGs. 6(b) and 6(d) show the spectrum of this system of spins under the 
influence of the large polarizing magnetic field and the field from the ferromagnetic 
sphere for two imaging angles, 0° and 45°, respectively. There are several 
observations with respect to FIGs. 6(a)-(d) that are pertinent to the image 
5 reconstruction process. It is immediately apparent that the projection along a single 
angle is not an ideal projection as described by Equation (3), but is modified by two 
distorting factors. The first distorting factor is that the projection in tomographic 
magnetic resonance microscopy of a single spin is a convolution of the point 
projection with the Lorentzian line shape L(q) along the projection line 

10 R,(g)= f . JJfoTIfo-fW [Equation (9)] 

[0062] The second, and more serious distortion, is the nonlinearity of the projection 
due to the nonlinear variation of the magnetic field along the z axis. As described 
above, in the configuration shown in FIG. 4, the one-dimensional radon transforms 
are not obtained by the uniformly separated parallel rays of constant Bz, but by the 
15 contours that have an increasing separation along the z direction, as FIG. 4(c) 
illustrates. 

[0063] FIG. 7(a) shows the projection at 45° superimposed by the ideal projection 
that one would wish to observe where the response of a spin is a delta function. As 
opposed to detecting a series of equally spaced delta functions whose heights 
20 represent the number of spins along a projection line, the image is a series of 

Lorentzian line shapes that are not evenly spaced and therefore do not represent a true 
sample projection. However, despite this distortion, if image reconstruction is 
performed by the Fourier transform filtered backprojection algorithm for parallel 
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projections of Equation (4), one obtains a reasonable reconstruction of the sample 
spin density, as FIG. 7(b) illustrates. The image in FIG. 7(b) was reconstructed from 
90 simulated projections at 2° apart without linearization and deconvolution. The 
image is mostly distorted on the edges of the sample where the nonlinearity most 
5 significantly effects the projection. In addition, since the projection of a single spin is 
not an ideal delta-function projection, but is a projection convoluted by the Lorentzian 
line shape (Equation (9)), the contrast in the image is also reduced as compared to the 
original phantom. 

[0064] The magnetic field component (Equation (5)) used for magnetic resonance 
10 imaging may be used to correct the nonlinear projection distortion prior to the initial 
step in the Fourier transform filtered backprojection algorithm. The resulting 
linearized representative projection at an angle of 45° and the full reconstruction from 
90 simulated projections at 2° apart are shown in FIGs. 8(a) and 8(b), respectively. 
As illustrated, the computed image is now a very good reconstruction of the original 
1 5 phantom. The reconstructed image still has worse contrast than the original due to the 
Lorentzian line shape convolution effect. However, if the single spin response is 
known (which may be assumed for the Lorentzian response for the single spin 
magnetic resonance) it is still possible to further improve the image reconstruction by: 
(a) first deconvolving each projection from the Lorentzian line shape, then (b) 
20 perform the linearization, and finally (c) perform the Fourier transform filtered 
backprojection by using the algorithm for parallel projections of Equation (4) to 
obtain the final image. 

[0065] A deconvolved and linearized projection at a representative angle of 45° is 
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shown in FIG. 9(a), and the final computed image from 90 simulated projections at 2° 
apart is shown in FIG. 9(b) where excellent reconstruction of the original image is 
observed. 

[0066] As described above, in accordance with one or more embodiments of the 
5 invention, two-dimensional tomographic imaging can be performed on samples that 
are different from the regular ordered arrays of spins. FIGs. 10(a)- 10(c) illustrate a 
simulated image where the term "cal tech" appears. The projections of the sample 
will not have regular periodic features as presented above. Nevertheless, the 
presented tomographic magnetic resonance imaging method is designed to 

10 successfully reconstruct the original image, provided a sufficient number of 
projections are obtained. FIGs 10(a)-10(c) show the computed images of three 
successively improved reconstructions, respectively, where 90 simulated projections 
at 2° apart are used. FIG. 10(a) is a direct reconstruction using Equation (4) without 
performing the deconvolution and linearization of projections. As previously 

15 described, the image has significant distortions on the edges due to the non-linear 
distortion of projections. 

[0067] FIG. 10(b) is a reconstruction of linearized projections that have not been 
deconvolved from the Lorentzian line shape, and improved reconstruction is 
observed. FIG. 10(c) is a reconstruction that uses both deconvolution and 
20 linearization of projections before the final computed reconstruction using Equation 
(4). In FIG. 10(c), excellent reconstruction of the image as well as the contrast are 
illustrated. 

[0068] As described above, the sample is located at the position where the contours 
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of constants are perpendicular to the sample surface, and sample placement with 
respect to the ferromagnetic sphere. FIG. 1 1 illustrates several other possible sample- 
sphere orientations and the respective sample rotation axis required for computerized 
tomographic image reconstruction. In all four cases shown in FIG. 1 1, the contours of 
5 constant Bz are perpendicular to the sample, but vary slightly in the intersecting radius 
of curvature. Which relative sample-sphere orientation is chosen will depend on the 
specific experimental and space requirements that might depend on the ferromagnetic 
probe [45][79], detector (e.g., mechanical force detectors [80]), rf coils [81], and 
polarizing magnet configurations. Detailed analysis of the effect of different 
10 configurations on the resolution, sensitivity, linearity, and sample-detector coupling, 
as well as the potential of extending the technique to three-dimensional computerized 
tomographic imaging may also be conducted. 

[0069] In view of the above, it may also be observed that in most room temperature 
nuclear magnetic resonance (NMR) experiments, the fractional polarization is quite 
15 small, on the order of 10" 5 , and improves at helium temperatures to values on the 

-3 

order of 10" . Such fractional polarization may have important implications for 
observing the magnetic resonance signal at those temperatures due to the significant 
spin noise background. However, cantilever operation at millikelvin temperatures has 
been demonstrated [18] and similar technical advancements might alleviate the spin 
20 noise problem. Nonetheless, it should be noted that the spin noise problem in 
magnetic resonance is a topic of ongoing investigations and still somewhat 
controversial and not well understood [47] [48]. 

[0070] Additionally, the invention has focused on cobalt as the material of choice, 
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but magnetic materials with similar saturation magnetization but higher anisotropy, 
such as rare earth alloy PrFeB, may be used to reduce thermal fluctuations of the 
sphere probe [20] [49]. 

[0071] It should also be emphasized that although the invention may be aimed at the 
5 MRFM community for the realization of tomographic magnetic resonance imaging, 
many other detection systems (e.g., to improve imaging capabilities in the detection of 
a small number of spins) may also use the invention. Examples of such alternative 
systems may include the detection of direct transfer of angular momentum [82][83] 
and energy [84] [85] to the spin population in magnetic resonance using 

10 micromechanical cantilevers, flux-detection class of magnetic resonance detection 
schemes that may include microcoil NMR [50] [87], microsuperconducting quantum 
interference device detectors [51][88], Hall sensors [52], superconducting resonators 
[53], and microwave waveguides [54], which may provide additional routes to the 
proposed potentially atomic resolution magnetic resonance tomographic microscopy. 

15 Flux detectors might also have additional advantages when compared to the force 
detectors with regards to the linearity of the probe-sample coupling [89], and specific 
sample-detector coupling mechanisms [25] [89] [90] may also be used. 
[0072] Embodiments of the invention may not consider the specific phase 
coherence between the spins, as it is yet undetermined if the technique would be 

20 performed in cw [55] or pulsed nuclear magnetic resonance mode [56]. The presence 
of unusually large magnetic field gradients described above may have significant 
effects on the nuclear magnetic resonance pulse sequence used, but research into this 
field gradient regime may continue [57][58] and may result in the development of 
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pulse sequences further that will make this proposal feasible in accordance with the 
invention and the pulsed NMR mode. 

[0073] In addition to the above, one or more embodiments of the invention integrate 
a ferromagnetic sphere with an inductive micro-coil that may be used to detect and 
5 image spin distribution by the tomographic method described above. The history of 
custom wound micro-coils is well established, with examples ranging from the 
measurements of the Aharonov-Bohm effect [64], propagation of magnetic domain 
walls [65], high-sensitivity detection of electron [66][67], and nuclear spin resonance 
[68], microscopic electron and nuclear spin magnetic resonance imaging [69][70][71], 

10 magnetic micromanipulation [72] -[74], and the application of high gradient 

alternative magnetic fields in the high-sensitivity force magnetometry [75]. Further 
techniques fabricate planar or three-dimensional coil structures [76] -[78]. 
[0074] One or more embodiments of the invention utilize a micro-coil that is made 
by winding a 30 jim diameter magnet wire around an optical fiber. The functionality 

15 of such micro-coils has been extended towards magnetic resonance tomographic 

microscopy by incorporating a 15 jam diameter ferromagnetic Nickel micro-sphere on 
the cleaved end of the optical fiber. The coil is positioned near the ferromagnetic 
sphere on the cleaved end of the optical fiber, as FIG. 12 shows, to maximize the 
detection of the time varying magnetic field from the region of the sample in magnetic 

20 resonance. In FIG. 12, a nickel ferromagnetic sphere, 15-/im in diameter, is mounted 
on the cleaved end of a 65-/xm diameter optical fiber. A magnet wire, 30-/Ltm in 
diameter, is wound around the fiber for RF transmission and magnetic resonance 
detection. 




Potential Devices 

[0075] The invention may be implemented in several potential devices including 
inductive micro-coil structures, as well as nano-mechanical resonators that utilize 
5 nano-wire resonators integrated with colloidal metal nano-reflectors for optical 
detection of vibration. 

[0076] For the inductive detection of magnetic resonance, two devices may be used. 
One involves integrating ferromagnetic spheres with micro-coils wound on optical 
fibers (illustrated in FIG. 12), while the other one utilizes micro-coil structures wound 

10 on the capillary tubes filled with ferromagnetic material (FIG. 13). 

[0077] For the case of mechanical detection of magnetic resonance vibration 
through optical means, nanowire resonators may contain all the necessary components 
for force detection of magnetic resonance: a) nano-magnetic particle for applying 
ultra-high magnetic field gradients, b) plasmonic metal colloids for optical detection 

15 of resonator vibration, and c) non-magnetic, non-plasmonic nano-wire structure as a 
cantilever for high sensitivity detection. FIGS. 14(a) and (c) illustrate electron 
microscope images of a multi-sectioned nanowire resonator that may be used for high 
sensitivity detection in accordance with one or more embodiments of the invention. 
In this regard, FIG. 14(c) illustrates a magnetic structure (on the left part of the 

20 image), a silver ball, and a non magnetic structure (on the right part of the image). 
FIG. 14(b) illustrates an optical microscopic image of the device illustrated in FIG. 
14(a). FIG. 14(b) illustrates how bright the ball looks in an optical image. FIG. 14(d) 
illustrates the method of suspending the nanowire resonantor of FIGS. 14(a), (b), and 
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(c). As illustrated, silicon is grown on a glass slide followed by a nanoresonator 
deposition. Electron beam lithography is then used to define clamping pad on the 
nanoresonator. Metal is desposited and lifted off followed by silicon etching and 
nanowire suspension. 

5 [0078] FIG. 15 illustrates an interferometric optical detection set-up (containing the 
nanowire resonantor fabricated in FIGS. 14(a)-(d)) that could be used for the 
measurement of a nano-mechanical system. The nano-resonator 1502 is suspended 
between glass slide surfaces 1504 within a vacuum sealed chamber 1506 by the O- 
rings 1508 and vacuum valves 1510. An incoming laser beam 1512 is split by a beam 
1 0 splitter 1514. The reference signal 1516 is/may not be used and the data signal 1518 
is focused on the colloidal nano-reflector (of the nano-resonator 1502). The data 
signal 1518 is focused by using optical lens 1520 and immersion oil/plasmon resonant 
colloid 1522 (that prevents a reflection from the bottom surface of the bottom slide 
1504). 

1 5 [0079] The reflected beam from the plasmon resonant colloid and the reflected 

beam from the glass slide surface form the interferometric signal 1524 that is used for 
the measurement of the nano-mechanical resonator 1502 vibration. One can also 
actuate the resonator by magnetic means using the micro-coil fields/microcoil actuator 
1526. Thus, FIG. 15 presents a device for detecting movement through optical means 

20 in accordance with one or more embodiments of the invention. 



Logical Flow 

[0080] FIG. 16 is a flow chart illustrating the logical flow for conducting magnetic 
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resonance tomographic microscopy in accordance with one or more embodiments of 
the invention. At step 1602, a non-crystalline sample (that may be two-dimensional 
or three-dimensional) is placed under the influence of a static polarizing magnetic 
field. Such a sample may take many forms including a small protein or molecule, cell 
5 or other biological or non-biological material of interest. At step 1604, a radio- 
frequency field is introduced perpendicular to the static polarizing first magnetic field. 
[0081] At steps 1606 and 1608, two or more magnetically resonant spins of the 
sample are simultaneously obtained. As part of this process, the sample is placed 
under the influence of a ferromagnetic sphere having a second magnetic field at step 
10 1606. Such a ferromagnetic sphere may be -10 times the size of the sample. The 
placement of the sample such that contours of the second magnetic field from the 
ferromagnetic sphere are perpendicular to the sample in the x-z plane, but varying 
slightly in an intersecting radius of curvature. 

[0082] At step 1608, the sample is sequentially angularly rotated, around a 
15 prescribed axis, with respect to the ferromagnetic sphere. Thus, multiple image 
projections of spins are obtained by rotating the sample with respect to the sphere. 
[0083] At step 1610, an image of the sample is reconstructed based on the imaged 
spins projections using computerized tomography. Such reconstruction may be 
performed using a Fourier transform filtered backprojection algorithm for parallel 
20 projections of the magnetic resonance spins from the ferromagnetic sphere through 
the sample. In this regard, the reconstructing may be performed by a Fourier 
transform of a radon transform, multiplication by a ramp function in conjugate space, 
followed by an inverse transformation, and integration over all angles in accordance 
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with Equation [4] above. 

[0084] To more clearly reconstruct the image, the invention may correct for 
nonlinear projection distortion prior to a first step in the Fourier transform filtered 
backprojection algorithm. Alternatively (or in addition), each projection may be 
5 deconvolved from a Lorentzian line shape and linearized prior to performing the 
Fourier transform filtered backprojection. 

[0085] Thus, the above described methodology provides computerized tomographic 
microscopy in magnetic resonance imaging using the magnetic fields and field 
gradients from a ferromagnetic probe. A configuration is utilized where a two- 

10 dimensional sample is under the influence of a large static polarizing field, a small 
perpendicular radio-frequency field, and a magnetic field from a ferromagnetic 
sphere. Despite the non-uniform and non-linear nature of the fields from a 
microscopic magnetic sphere, the concepts of computerized tomography can be 
applied to obtain proper image reconstruction from the original spectral data by 

15 sequentially varying the relative sample-sphere angular orientation as illustrated in 
FIG. 11. 

[0086] In the invention, atomic resolution magnetic resonance imaging of discrete 
periodic crystal lattice planes using ferromagnetic probes are extended to two- 
dimensional imaging of non-crystalline samples with resolution ranging from 
20 micrometers to Angstrom scales. 
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Conclusion 

[0176] This concludes the description of the preferred embodiment of the invention. 
One or more embodiments of the invention provide a technique of two-dimensional 
10 tomographic magnetic resonance imaging using ferromagnetic probes that is distinctly 
different from other scanning probe techniques that achieve surface atomic resolution 
imaging, such as scanning tunneling microscopy [59], atomic force microscopy [60], 
spin polarized tunneling microscopy [61], and spin resonance tunneling microscopy 
[62][63]. 

15 [0177] In prior art methods, the images are obtained by raster scanning an 

atomically sharp probe over a two-dimensional surface of the sample, and at each 
pixel of the image a measured parameter such as a tunneling current or atomic force is 
displayed. In the two-dimensional magnetic resonance tomographic reconstruction 
microscopy of the invention, the ferromagnetic probe has no physical atomic scale 

20 features to directly achieve atomic resolution, but rather provides the magnetic field 
gradients that are sufficiently large to atomically resolve magnetic resonance of 
individual spins along only one direction. 

[0178] Embodiments of the invention do not rely on a point-by-point data 
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acquisition for two-dimensional imaging, but rather is conducted by sequential 
angular rotation of the sample with respect to the ferromagnetic probe and external 
field around a prescribed axis. At a single angular orientation and magnetic field 
value many magnetically resonant spins of the sample are detected simultaneously. 
However, by recording such a multispin signal at many angular orientations and 
magnetic field values, it is possible to reconstruct the two-dimensional image 
representation using the methods of computerized tomography. 
[0179] The foregoing description of one or more embodiments of the invention has 
been presented for the purposes of illustration and description. It is not intended to be 
exhaustive or to limit the invention to the precise form disclosed. Many 
modifications and variations are possible in light of the above teaching. It is intended 
that the scope of the invention be limited not by this detailed description, but rather by 
the claims appended hereto. 
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